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Abstract: Our study evaluated the effects of environmental variables on the assemblages of the sub- 
order Zygoptera, and tested the hypothesis that environmental variables are more important determi- 
nants of the structure of these assemblages than limnological variables in streams. We sampled 17 
streams in the Carajás National Forest and tested our hypothesis using a linear regression analysis, 
with the zygopteran species composition, richness, and abundance as the response variables. Our 
findings indicate that both limnological and physical variables influence, independently, the charac- 
teristics of the zygopteran assemblages. The riparian forest maintains the stability of the environ- 
ment and provides dispersal corridors, along which the zygopterans can reach alternative, suitable 
environments. The small scale of this study also implies that the continuity of the vegetation is es- 
sential for the dispersal of the zygopterans among different landscapes. The high levels of abundance 
recorded in the better-preserved environments may reflect the maintenance of specific habitats and 
resource availability. Riparian forest is crucial to the ecological equilibrium of the stream systems, 
although further research at a broader spatial scale that focuses on a greater diversity of variables 
should provide more robust insights into the phenomenon. 
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Introduction 


The distribution of aquatic insects has often been explained in terms of the influence of the physical 
structure of the habitat and the chemical composition of the water (Heino et al., 2005; Landeiro et al., 
2011; Landeiro et al., 2012; Cunha et al., 2015). In their seminal work, McArthur & McArthur (1961) 
proposed that habitat heterogeneity and its physical components have an important influence on the 
characteristics of the community, providing the physical conditions and resources required by the 
organisms to survive in the environment. The structural components of stream systems, such as 
woody and leafy debris, shape the distribution of their aquatic fauna (1.e., fish and benthonic inverte- 
brates), creating a heterogeneous gradient that determines the variation in the breeding and foraging 
patterns of different organisms (MacArthur & MacArthur, 1961; Thomaz & Cunha, 2010; Kovalenko 
et al., 2012). Environmental factors, such as ambient temperatures, the chemical compounds dis- 
solved in the water, and its turbidity, may also be important, functioning as environmental filters, to 
determine which species inhabit a given environment (Oliveira Jünior et al., 2015; Brasil et al., 2018). 


In the Amazon rainforest, the principal impacts on stream systems are derived from human activities, 
in particular, cash cropping (Cunha et al., 2015), cattle ranching, logging (Calvào et al., 2016), hydro- 
electric dams, and mining (Monteiro, 2005; Fearnside, 2006; Luiza-Andrade et al., 2020), all of which 
may influence the physical structure of the stream margins and the quality of the water, which affect 
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the stability of the environmental (Nessimian et al., 2008; Brasil et al., 2014a; Cunha et al., 2015; 
Oliveira Junior et al., 2015). Mining is widespread in the Amazon region, and can have a significant 
impact on many of the natural characteristics of the stream systems (Callisto et al., 1998). Mining may 
also affect specific features of the water, such as its pH and the concentrations of some chemical com- 
pounds, and any change in the input of energy may have unpredictable consequences for the resident 
fauna (Callisto et al., 1998; Enriquez-Espinosa et al., 2020). In Brazil, mining operations are covered 
by specific environmental legislation (federal law 12,561/2012), which legislates the preservation of the 
riparian vegetation of the streams within the area affected by the activity (Soares-Filho et al., 2014). 


One of the principal impacts on stream systems is the removal of the riparian vegetation (Lampert & 
Sommer, 2007; Nessimian et al., 2008; Brasil et al., 2014b). Nessimian et al. (2008) proposed the 
Habitat Integrity Index (HII), which quantifies the stability and integrity of this vegetation, and pro- 
vides a systematic estimate of its physical structure. The riparian vegetation plays a vital role as a 
saturation zone, by controlling the input of sediments into the stream bottom (Callisto et al., 2001; 
Astudillo et al., 2016), and maintaining the equilibrium of the environment by regulating air and 
water temperatures, the incidence of sunlight, shading, and the input of the allochthonous materials 
that support the system's trophic webs (Maloney & Weller, 2011; Oliveira Júnior et al., 2015; Mon- 
teiro-Junior et al., 2015; Carvalho et al., 2018). 


Environmental disequilibrium caused by the suppression of the riparian vegetation may have a major 
impact on the fauna of aquatic systems (Nessimian et al., 2008). A number of ecological studies have 
demonstrated that these impacts harm the natural assemblages found in the aquatic systems (Brasil 
et al., 2014a; Oliveira Junior et al., 2015; Siquiera et al., 2015; Astudillo et al., 2016; Carvalho et al., 
2018), in particular, by decreasing both species diversity and abundance. The removal of the riparian 
vegetation impacts environmental stability primarily by altering ambient temperatures, water qual- 
ity, and the shade provided by the vegetation (Carvalho et al., 2013). Impacts on the riparian vegeta- 
tion also modify the connectivity between streams, which alters dispersal patterns and the distribu- 
tion of organisms at both local and regional levels (Carvalho et al., 2013; Heino, 2013; Monteiro 
Junior et al., 2015; Tonkin et al., 2017). 


The suborder Zygoptera (Odonata; Insecta) comprises insects that are primarily specialists adapted 
for specific environmental conditions, typically within a narrow range (Corbet & May, 2008; Olivei- 
ra Junior et al., 2015). Corbet (1999) classified most zygopteran adults as thermal conformers, with a 
high thermal conductance, which means that their body temperature is determined by their immedi- 
ate environment. Thermal conformers may be affected profoundly by modifications in the riparian 
forest, and are thus considered to be habitat specialists, which supports their role as bioindicators of 
environmental impact (Oliveira Junior et al., 2015; Calvão et al., 2016). In most cases, zygopteran 
species richness and abundance are reduced, and in extreme cases, local extinctions may occur 
(Oliveira Jünior et al., 2015; Brasil et al., 2018; Mafuwe & Moyo, 2020). Even when the riparian forest 
is not affected directly by the anthropogenic impact, the conditions of the surrounding matrix may 
also be important for the dispersal of the insects beyond the riparian zone (Calvao et al., 2016), which 
creates a habitat gradient that the species need to traverse (Dolny et al., 2014). In this scenario, 
modifications in ambient temperatures or the vegetation cover over the streams would have negative 
impacts on adult zygopterans, which require very specific conditions (Corbet, 1999; Monteiro Jünior 
et al., 2015). Understanding the influence of environmental factors and connectivity on the dynamics 
of stream assemblages is a fundamental objective in aquatic ecology. Adult zygopterans are good 
models for the assessment of anthropogenic impacts (Oliveira Junior & Juen, 2019) through the as- 
sessment of dispersal patterns in relation to specific environmental interruptions (Carvalho et al., 
2013; Juen et al., 2014; Oliveira Junior et al., 2015). 


In the present study, we evaluated the role of physical variables, based on the HII and the features of 
the nearest mining site (pit size and distance), and limnological variables on the characteristics of the 
zygopteran assemblages observed at the sampling sites, located in stream systems in the Carajas 
National Forest, in northern Brazil. We hypothesized that elements of the physical structure structure 
of the streams are more effective drivers of zygopteran diversity (species composition, richness, and 
abundance) than limnological variables, such as dissolved gases and solids. In this case, the HII pro- 
vides a surrogate of physical structure, and higher HII scores indicate better preserved sites with a 
high degree of environmental equilibrium, which is essential for the maintenance of diversity in both 
the streams and the surrounding habitats. In general, the integrity of streams will decline with in- 
creasing proximity to mining operations (Enriquez-Espinosa et al., 2020). 
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Materials and Methods 


Study area 


This study was conducted in the Carajás National Forest, in southeastern Pará state, northern Brazil 
(06°0’00” S, 50?24'0" W; Figure 1). This region has a humid tropical climate with dry winters, clas- 
sified as AWi in the Kóppen system (Peel et al., 2007). More than 95% of this conservation unit is 
covered by dense hardwood forests, with the remaining area consisting of a type of savanna, known 
locally as the “canga”, which grows on soils rich in iron-ore, and is referred to as the Iron Triangle 
(Ab'Saber, 1986; Luiza-Andrade et al., 2020). The rainy season lasts from November to May, and the 
dry season, from June to October, with total annual precipitation of 1800—2300 mm (Moraes et al., 
2005). The study area is 400—900 m above sea level (a.s.1.), and contrast with the surrounding area, 
which has altitudes of between 80 m and 300 m a.s.l. (Soares-Filho et al., 2014). 


South America State of Pará 


«7 


LIE € D, i 
2 E. 
f 


es 4 180 360 Km| 
——— 

P d 
"e 


^ Conservation Units 


is » 
oe d SSIES 
Ane Pe ia 
2 TSA or 
 TRPIIAPEAQUIAI FLONA wv 
s 
x TIACAIUNAS FLONA f MU 
Sole iw CWARAPEGEEADQAPA 


Apel. 


A 
Eee NE 


CAMPOS FERRUGINOSOS LN 
( = 
A 


0 125 25Km 
LLL 


[ER Forest a Mining 
eS Canga eS Hydrography 
qe Regeneration Cloud 4 Sampled sites 





Figure 1. The 17 study streams sampled during the present study in the Carajás National Forest, Para 
state, Brazil. The land use categories are shown in the legend. 


The principal local hydrographic basin includes the Itacaiinas and Parauapebas rivers, and covers an 
area of more than 11,500 km?, which includes several first-order tributaries (Soares-Filho et al., 2014). 
The streams in the Itacaiünas watershed have white (milky coffee-colored) water during the dry sea- 
son, changing to yellow or red in the rainy season, when the sediment load increases (Ruivo & Sales, 
1989). Iron ore has been extracted within this area since 1970, which prompted the creation of the Na- 
tional Forest (IBAMA, 2003; Luiza-Andrade et al., 2020). As a result, the companies operating within 
this area must strictly respect Brazilian environmental legislation (the forest code, federal law 
12,561/2012), which mandates the preservation of a 30-m buffer of riparian vegetation around all 
streams (ITV, 2017). 
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Biological sampling 


Samples were collected from 17 study streams in October and November 2015, which corresponds to 
the transition from the rainy to the dry season. The adults were collected between 10:00 am and 2:00 
pm, the period of the greatest incidence of sunlight, which is the ideal condition for the physiological 
requirements of the zygopterans (Corbet, 1999; Oliveira Junior et al., 2015). Adult zygopterans were 
collected at each study stream within a standard 100 m stretch divided into 20 sections of 5 m. Within 
each section, 10 transects were demarcated and named from A (downstream) to K (upstream). The in- 
sects were collected by fixed-area sweeping (Juen & De Marco, 2011; Pinto et al., 2012; Monteiro 
Junior et al., 2015) with an entomological net (diameter: 40 cm; depth: 65 cm) attached to a 90-cm 
aluminum pole. Once collected, the specimens were stored in paper envelopes and immersed in ace- 
tone P.A. (pure for analysis) for 12 hours, following the protocol of Lencioni (2006). The specimens 
were identified at the Ecology and Conservation Laboratory of the Federal University of Para in Belém 
using the relevant taxonomic keys (Garrison et al., 2006, 2010; Lencioni, 2005, 2006). The biological 
data were collected during the field campaigns of the integrated program for the monitoring of bioindi- 
cators in the Carajas National Forest, which aims to mitigate the impacts of the local mining operations 
on the region’s biodiversity. Fieldwork at the Carajas Complex of the Vale S.A. mining corporation was 
authorized by the Brazilian Institute for the Environment and Renewable Natural Resources (IBAMA). 


Environmental variables 


The environmental dataset was composed of 12 variables (Supplementary Material; Tab. 2), which are 
all important drivers of zygopteran diversity, by determining the capacity of these insects to establish a 
viable population (Assis et al., 2004; Oliveira Junior et al., 2015, Calvào et al., 2016; Oliveira-Junior & 
Juen 2019). The HII (Nessimian et al., 2008) assesses the physical structure of the stream and its sur- 
roundings, and was calculated based on the values of the 12 environmental variables recorded at each 
stream. The scores for each variable in this protocol vary from 0 to 1, with the values closer to 0 repre- 
senting more impacted environments, and those closer to 1 representing more pristine habitats. A num- 
ber of previous studies have applied this index to the analysis of the variables that influence the occur- 
rence and distribution of odonates, and have provided important insights into the diversity of the 
zygopterans (Monteiro Jünior et al., 2014; Oliveira Jünior et al., 2015). In a recent meta-analysis, Brasil 
et al. (2020) obtained excellent results using the HII to evaluate the structure of communities of aquatic 
insects, and confirmed that this index is an effective parameter for the assessment of the environmental 
impacts that affect the aquatic biota. We also included the Impact Index (II) in our analysis, which was 
provided by the company that manages the study area. This index (II) measures the reduction of the 
habitat, considering the proximity of the mining activities to the sampling sites, with the extent of the 
impact on the fauna being expected to increase with increasing impact. On the other hand, the impact 
on the streams will be inversely related to the squared distance from the nearest mining zone (Espinoza 
et al., 2020). In the present study, we refer to the physical-chemical parameters as the limnological vari- 
ables and the HII, II, and the distance to the nearest mining pit as the physical structure variables. 


Data analysis 


We applied a Variance Inflation Factor (VIF) to the limnological variables (James et al., 2014) to 
control for multicollinearity, and omitted the metrics with values that exceeded 5. This left three 
variables (Dissolved Oxygen = DO; environmental temperature = ENVT; turbidity — turb), which 
were included in the model, together with the HII and II (Supplementary Material; Tab. 2). A Princi- 
pal Components Analysis (PCA) was applied to minimize the dimensionality of environmental data 
and describe the distribution of the sampling sites, for which the variables measured on different 
scales were standardized (Legendre & Legendre 2012). The broken stick criterion was applied to the 
data to identify the most important axes derived from the variables, considering only the values that 
were at least as large as those predicted by the model (Jackson 1993). A Principal Coordinates Anal- 
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ysis (PCoA) was applied, based on the Bray-Curtis dissimilarity index, to verify the variation in 
species composition among the study streams (Legendre & Legendre 2012). For this, species abun- 
dance was log (Log [x*1]) transformed, and the first two axes were retained (Legendre & Legendre 
2012). To test our hypothesis, the first PCoA axis was used as a surrogate of species composition. 


Limnology + 
Physical Zygoptera 
data 
elements ae 
VIF (=<5) SFA COSI 
Bray-Curtis 
i + 

PCA Remained STA A 


variables PCOA cnaig 


Best subset regression waco 


(Zygoptera ~ Limnology /physical elements) 


Figure 2. A summary of the statistical procedures used in the present study. AAICc = difference in the 
Akaike Information Criterion in comparison with the lowest AlCc value; HIl= Habitat Integrity Index; PCA 
= Principal Components Analysis; PCoA = Principal Coordinates Analysis; N = zygopteran abundance; S = 
zygopteran species richness. 


We ran a Multiple Linear Regression (MLR) to assess the responses of the zygopteran diversity met- 
rics (species richness, composition, and abundance) to the predictor variables (limnological param- 
eters, physical structure (HII and ID), and the distance from the streams to the mining pits) using best 
subsets regression (Harrell 2001). This procedure involves compiling the regression models of all the 
possible combinations of predictor variables and selecting the best model, which meets the pre-de- 
fined criteria (Lynch et al., 2018). The models were ranked using the Akaike Information Criteria for 
small sample size (AICc), to identify the optimum models (Burnham & Anderson 2002). The best 
model from each dataset had the lowest AICc value, while models with AICc values of less than 2 
(AAICc<2) were also considered to be good predictors of diversity. The PCA and PCoA were run in 
the R environment (R Core Team 2017), using the ‘vegan’, ‘MASS’, and ‘ggplot2’ packages (Wick- 
ham 2016), and the best subsets regression procedure was run in Systat 13. The analytical procedures 
are summarized in Figure 2. 


Results 


General description of the zygopteran assemblages 


A total of 343 zygopteran adults were collected during the present study, representing four families 
(Coenagrionidae, Polythoridae, Calopterygidae, and Megapodagrionidae), and 19 species (Supple- 
mentary Material; Tab. 1). The most abundant species was Chalcopteryx rutilans Rambur 1842, with 
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104 individuals. The second most abundant species was Argia tinctipennis Selys 1865, with 69 indi- 
viduals, followed by Hetaerina indeprensa Garrison 1990 with 35 individuals, and Mnesarete aenea 
Selys 1853 with 23. 


General description of abiotic variables and species composition 


The PCoA explained 46.3% of the total variance on its first two axes, with 28.53% in the first axis, 
and 17.77% in the second (Figure 3A). However, the analysis did not the existence of a defined distri- 
bution pattern in the species composition, which varied randomly. 


The broken stick model retained both of the first two PCA axes, which together explained 69.58% 
of the total variation, with the sites with the highest II scores appearing the lower left corner of the 
ordination, and some of the sites with higher HII scores in the upper left corner (Table 1; Fig- 
ure 3B). Dissolved oxygen was the variable that contributed most to the first axis (0.747), while 
turbidity (—0.919) and the II (—0.561) had the greatest negative influence on the first axis (Figure 
3B). The HII was the most important variable on the second axis, with a positive score of 0.899 
(Figure 3B). 
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Figure 3. Ordination of the abiotic and biotic data - A: Principal Coordinates Analysis (PCoA) applied to 
the data on the zygopteran species composition; B: Principal Components Analysis (PCA) applied on the 
abiotic data. Turb = turbidity; DO = dissolved oxygen. HII = habitat integrity index; Il = impact index. 
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Table 1. Mean, standard deviation (SD), and maximum/minimum values of the five variables included in 
the PCA (Principal Component Analysis) of environmental factors. The values were obtained at the study 
streams in the Carajas National Forest, Para state, Brazil. 




















Variables Mean SD Maximum Minimum 
Dissolved Oxygen (mg/L) 7.64 1.27 8.60 3.48 
Environmental Temperature, 27.48 6.72 30.70 25.00 
ENVT (°C) 

Turbidity (NTU) 5.47 4.24 13.90 0.67 
Habitat Integrity Index (HII) 0.87 0.11 0.56 0.97 
Impact Index (II) 10.81 15.86 60.04 0.13 




















Table 2. Multiple linear regressions (best subsets) of the zygopteran diversity metrics and their relation- 
ship with the environmental and physical predictors at the study streams in the Carajás National Forest, 
Para, Brazil. AIC = Akaike Information Criterion; AlCc = Akaike Information Criterion corrected for small 
samples; AAICc = AICc ranking based on the minimum value; DO = Dissolved Oxygen; ENVT = environ- 
mental temperature ; HII = Habitat Integrity Index; PCoA1 = first axis of the Principal Coordinates Analy- 
sis; TURB = turbidity; f = beta coefficient. 















































Selected model AIC AICc AAICc | Adjusted R? | p 
Species richness 
DO + ENVT 77.163 80.496 | 0 0.231 —0.834; 0.592 
DO + TURB 78.131 81.464 | 0.968 0.186 —].191; -0.279 
HII 79.622 81.468 | 0.972 0.067 7.195 
ENVT 79.694 81.540 | 1.044 0.063 0.435 
Abundance 
HII 131.238 | 135.084 | 0 0.099 39.42 
DISTANCE 133.980 | 135.827 | 0.743 0.058 9.389 
DISTANCE + TURB 132.750 | 136.083 | 0.999 0.166 0.002; 1,238 
ENVT * TURB 133.025 | 136.359 | 1.275 0.152 4.409; 1.701 
Composition (PCoA1) 
DO 4.078 5.924 | 0 0.037 —0.063 
DISTANCE 4.704 6.550 | 0.626 0.001 >0.001 
DO + TURB 4.374 7.707 | 1.783 0.066 —0.108; —0.021 
DO +DISTANCE 4.430 7.763 | 1.839 0.063 —0.069; 70.001 
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Effects of environmental variables and physical structuring on the 
Zygoptera 


The best subsets regression indicated that both types of environmental parameter influence all the 
zygopteran diversity metrics in some ways, although in most cases, the relationship was relatively 
discreet (Table 2). However, our hypothesis was supported, in part, by the results, given that the HII 
influenced only zygopteran abundance when analyzed on its own, while the II did not influence the 
diversity or abundance of the zygopterans (Table 2). Species richness was best explained by the 
model that combined DO and environmental temperature (R? adj = 0.231). However, the model that 
included only the HII best explained the abundance of individuals, albeit with a greatly reduced ex- 
planation value (R? adj = 0.099). The distance from the stream to the near mining pit was also preva- 
lent in the zygopteran abundance models, although once again, with only low values (Table 2). The 
best model selected to explain the zygopteran species composition was formed exclusively by DO, 
but once again, with a greatly-reduced value of explanation (R? adj = 0.037), with a similar pattern 
being observed in the distance from the pit (Table 2). It is important to note that turbidity was the 
variable selected most frequently by the AAICc (« 2), playing a role in the models of all the diversity 
metrics, followed by all the other variables, except the II (Table 2). 


Discussion 


The hypothesis tested in the present study was partially corroborated, by showing that the HII was 
the principal parameter of physical structure that drove zygopteran abundance. Although the limno- 
logical features had a considerable influence on all the zygopteran diversity metrics, and were the 
principal drivers of species composition and abundance, the low explanatory values do not uphold 
our hypothesis. The HII values, which represent the physical structure of the stream, also varied 
discreetly, from 0.57 to 0.99, which indicates that most of the study sites had a relatively high level of 
habitat integrity, with good vegetation cover and stable environmental conditions. This may have 
accounted for the broadly similar patterns found among sites, given that the relative homogeneity of 
the riparian vegetation is an important factor in zygopteran biology (Calvao et al., 2016; Carvalho et 
al., 2018). 


The limnological parameters provided important insights in the present study, given their fundamen- 
tal influence on the zygopteran communities (Mendes et al., 2018). Turbidity can influence the selec- 
tion of oviposition sites by the female, for example (e.g. Carvalho & Nessimian, 1999), given that 
they typically prefer more transparent water. The influence of the physical-chemical parameters of 
the water on the occurrence of adult zygopterans should be interpreted with caution, however, given 
that the adult and nymph stages are closely related (Valente-Neto et al., 2015; Mendes et al., 2017; 
Iwai et al., 2017), which implies that the factors affecting one stage may affect the other by default 
(Brito et al., 2020). As the present study focused on the adults, then, a more integrated approach 
would be needed to provide more conclusive insights into these patterns. Dissolved oxygen was the 
principal variable influencing zygopteran species richness and composition, with negative beta coef- 
ficient values indicating a certain loss of integrity (Calvào et al., 2018). This outcome is quite interest- 
ing because, even though the HII was not retained as the principal variable in all the models, the DO 
does show that changes in water composition are occurring and may influence the odonate larval 
stages (Mendes et al., 2018). We cannot make further inferences, however, given that we did not 
evaluate this life stage in the present study. 


The HII evaluates the completeness and integrity of the riparian forest in the context of the potential 
impact of human activities (Nessimian et al., 2008). The integrity of this habitat is fundamental to the 
maintenance of stable temperatures and the filtering of the sediment input, as well as providing a 
dispersal corridor for many stream-dwelling species (Oliveira Junior et al., 2015; Calvão et al., 2018; 
Carvalho et al., 2018). In the present study, in fact, only one stream (ssl) had an HII lower than 0.80, 
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even though some sites were located in the vicinity of mining zones. This undoubtedly reflects the 
systematic application of federal law 12,561/2012, which mandates the preservation of the riparian 
zone by companies extracting natural resources, in order to preserve ecological networks and protect 
water resources. These protective measures ensure the conservation of the riparian zone, an impor- 
tant factor influencing the diversity of zygopteran species (Monteiro-Junior et al., 2015; Nagy et al., 
2019), given that these organisms have highly specific physiological requirements, which demand not 
only vegetation cover, but also environmental stability. 


The present study focused only on zygopteran ectotherms, which are unable to regulate their body 
temperature (Oliveira Junior et al., 2015), which further reinforces the importance of the riparian 
influence zone for the occurrence of species. In addition, the distance from each stream to the the 
nearest mining pit was a factor in both zygopteran abundance and species composition, albeit in a 
secondary role, which indicates that this impact is not extensive enough to impede the dispersal of 
the damselflies across the environmental gradient. Despite being much smaller than the dragonflies 
(suborder Anisoptera), the damselflies appear to be capable of traversing the habitat gradient of the 
study area to reach other suitable sites. It is important to note here that the study area is located 
within a large broad conservation unit, and that even the streams in closet proximity to mining pits 
have intact riparian forest, which ensures the connectivity among and guarantees species dispersal. 
Alternatively, it is possible that the stability of the physical structure of the habitats ensures the ho- 
mogeneity of the chemical composition of the water, its dissolved oxygen concentration, and turbid- 
ity, which is essential for the zygopteran adults (Raebel et al., 2012; Calvào et al., 2016; Mendes et al., 
2018). 


Brasil et al. (2020) found that the HII is an effective predictor of the composition of aquatic insect 
communities in Brazilian stream systems. This study highlighted the correlation between environ- 
mental variables, such as dissolved oxygen and turbidity, and the HII, which quantifies the physical 
structure of the stream (Nessimian et al., 2008). Higher HII scores indicate a greater degree of integ- 
rity, that is, environmental stability. Habitat integrity thus has a fundamental influence on the envi- 
ronmental dynamics of these systems, including the chemical composition of the water and physical 
conditions such as air temperature, which are all crucial for zygopterans. The well-preserved riparian 
forest found throughout the area of the present study area may account for the abundance of some 
species, such as C. rutilans, whose females require specific breeding sites, which are only found in 
well-preserved environments (Resende & De Marco Jünior 2010). The C. rutilans females lay their 
eggs on fallen tree trunks and other structures derived from the riparian forest (Resende 2010), and 
their dependence on these structures would account for their marked association with more pristine 
environments (Resende & De Marco Jünior 2010; Calvào et al., 2016). As this species contributed a 
large proportion (3396) of the individuals collected in the present study, a greater distance from the 
zone of impact may have been necessary to guarantee the occurrence of this species, given that this 
distance was related to damselfly abundance. 


In the present study, both limnological factors (physical-chemical properties of the water) and 
physical structuring (HII and mining), had a significant influence on zygopteran diversity. That is, 
both factors mold the establishment, reproduction, and survival of these organisms. Environmen- 
tal stability is critical to the zygopterans, which are small insects with slender bodies. As for the 
habitat integrity index, the environmental features varied only discreetly among sites, indicating 
a high degree of stability, with low adjusted R? values. One other important factor here was the 
relatively reduced spatial scale of the study area, which likely reinforced the similarities among 
the sites. One way to increase the sensitivity of the analytical approach adopted here would be to 
compare a much broader selection of sites that encompass a greater degree of variation in both the 
HII and the environmental parameters. Overall, them the physical habitat appears to be the prin- 
cipal factor driving the structuring of zygopteran diversity, which reinforces the need to consider 
local factors for the planning of conservation measures for aquatic ecosystems in the Amazon 
region. 
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